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In ultrarelativistic heavy-ion collisions, in the 2 GeV< p± < Q GeV transverse mo- 
mentum region, the soft and semi-hard multiple scattering of the incoming nucleons in 
the nuclear medium results in the broadening of the expected hadronic (e.g. pion) p± 
spectra relative to proton-proton (pp) collisions. Thus, higher transverse-momentum re- 
gions are populated than in pp collisions. In a perturbative QCD based calculation we 
include the intrinsic transverse momentum {k±) of the partons in the nucleon (determined 
from pp collisions), augmented by the extra broadening obtained via a systematic analysis 
of proton- nucleus (pA) collisions in the energy range 17 < \fs < 39 AGeV. The original 
polynomial spectra are modified, and a nearly exponential spectrum appears in the region 
2^ P± ^ 3.5 GeV. At present RHIC energies (a/s =130 AGeV), the slope of the calcu- 
lated spectra is reminiscent of that of fluid-dynamical descriptions, but lacks any thermal 
origin. We determine and discuss the size of the modifications originating in multiple 
scattering, which lead to this state of affairs. 



1. Introduction 

With the start of RHIC operation new data are now available for high energy heavy 
ion collisions. This makes it possible to contrast the models to the experimental data and 
to learn about new effects in this so far unreachable domain. 

In this talk we discuss the features of some hadronic spectra in Au+Au collisions at 
130 GeV. First we give an introduction to the pQCD-improved parton model and fix 
our parameters in pp collisions. Next, we study the available pA experiments in order 
to understand the nuclear (Cronin) effect, the enhancement of particle production in a 
certain p±_ range compared to pp coUisons. Finally, we present our calculation for nucleus- 
nucleus (^4^4) collisions at present RHIC energy. 



2. Perturbative QCD 



The classical pQCD parton model assumes factorization of the hadron (meson M) 
production cross section, namely, the particle production is interpreted as a superposition 
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of elementary, partonic collisions 0, 

Em^^^ = / dXadXbdZc fa/p{Xa,Q^) fb/p{Xb, Q'^)'^^^^^ Dm/c{Zc,Q ) g^^^g_^ ^ ^^-^ 

where fa/p{xa, Q^) is the parton distribution function (PDF), i.e. the probability of finding 
parton a in the proton with momentum fraction Xa at transverse momentum scale Q, 
da/dt is the pQCD calculable cross section of quarks and gluons, while DM/cizcQ"^) is 
the fragmentation function (FF) of parton c into a meson M with the meson carrying a 
momentum fraction Zc of the parton at transverse momentum scale Q. The PDF and FF, 
contain all the non-perturbative information and are fixed by experimental data. In our 
calculation we use the leading order GRV PDF and the leading order KKP FF. 
(We also carried out approximate NLO calculations 0). 

Evaluating integral (^ numerically and comparing it to existing experimental data of 
charged and neutral pion and kaon production in pp and pp reactions in the range of 
\/s = 17 — 1800 GeV we found a factor of two discrepancy with the theory underes- 
timating the yield. Even more sophisticated NLO calculations failed to reproduce the 
pion spectra . This motivates us to introduce another non-perturbative parameter, the 
intrinsic transverse momentum of partons. Such a quantity is not fictitious: it has been 
measured experimentally 0. We consider the distribution of this intrinsic transverse mo- 
mentum to be Gaussian for simplicity, and fix its width from experiments. The integral 
for both incoming partons in Eq. (^ is modified to 

dXa fa/p{Xa,Q'^) > dXa d\±^a ^^^^ fa/p{Xa,Q^)- (2) 

The fit of the intrinsic transverse momentum width is shown in Fig. |l] for pions. A similar 
plot may be obtained for kaons (with much less experimental data available). Using the 
intrinsic transverse momentum distribution with a width given by the figure we were able 
to fit pp experimental data within 30% for p± > 2 GeV. 

Having fitted the pp data using Eqs. (0-|5) parameterized by the {k\)pp values, we now 
turn to the pA collisions. It was found in experiments that there is an extra nuclear 
enhancement in pA collisions compared to the simple scaling of pp data in a certain p_[_ 
range (Cronin effect) [§]. One needs to extend integral (|ll-@) as 



d p J dp 

Here tyi(6) is the nuclear thickness function at impact parameter b, tA^b) = J dzg{b, z) and 
the nucleon-nucleon cross section has a dependence on the width of intrinsic transverse 
momentum of the proton and nucleons in nuclei A, respectively. Furthermore, due to the 
nuclear environment the partonic distribution functions are modified with a shadowing 
and antishadowing region 

The theoretical interpretation of the Cronin effect is related to the nuclear environment, 
the average intrinsic transverse momentum of a parton increases due to multiscattering, 

{kl)pA = {kl)pp + C h{vA{b) - I), (4) 
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Figure 1. Average transverse momentum 
square fitted to pp experiments 



Figure 2. Slope parameter T, for scaled pp 
collision (dotted line) and central Au+Au 
collisions at 130 GeV (solid line) in vr pro- 
duction. 



where z/^(6) is the number of target nucleons in the channel swept by the incoming proton 
at impact parameter 6, h is an "effectivity" function, selecting the number of target 
nucleons contributing to the {k\) enhancement of the projectile and C is the average 
momentum square imparted in one nucleon-nucleon collision. 

For the effectivity function we tried a function linearly rising with the number of par- 
ticles inside the channel and saturating after n collisions. Requiring approximate target 
and energy independence of h we found that the best fit is achieved using n = 3 — 4 with 
average momentum square impart C„=4 0.4 GeV^. In the following we use this value 
for central nucleus-nucleus collisions. 

Next, we performed a simulation of the ongoing RHIC experiments on meson production 
at ^/s =130 GeV in Au+Au collisions. This requires a modification of Eq. (P) extending 
it to the case of two colliding nuclei, 

eJ< . /d^.dV UV)t,ii-r) ^„ WiUiiM. (5) 



(Tp J <rp 

Based on Fig. ^ we used {k\)pp = 1.6 GeV^. A simple calculation without the nuclear 
enhancement and shadowing agrees well with the peripheral PHENIX results on tt" pro- 
duction [|T^. However, for central collisions there is a disagreement between this pQCD 
model and the experimental data, indicating a strong suppression in the particle produc- 



tion. A possible candidate explaining this phenomenon is jet quenching |Tll . 

An interesting quantity to study is the slope of the meson production cross section. For 
peripheral collisions (no shadowing and nuclear enhancement) a pQCD calculation yields 
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a non-exponential spectrum with a slope parameter rising with the transverse momentum 
(Fig. dotted line). However, in case of central collisions in a p± window 2-3.5 GeV 
the slope parameter is almost constant resembling an exponential (thermal) spectrum. 
Predictions for lower transverse momentum are uncontrollable within pQCD. Our studies 
show that at SPS energies this plateau is more pronounced and extends to the highest 
transverse momenta measured. 



3. Conclusion 

In this talk we showed that at high transverse momenta {p± > 2 GeV) pQCD augmented 
with an intrinsic transverse momentum distribution of the partons gives a good description 
of pp, pA and AA collisions. A "thermal" slope is observed in central AA collisions in a 
transverse momentum window 2-3.5 GeV despite the non-thermal character of the model. 
At RHIC energies, however, a strong suppression is observed in hadron production in 
central collisions compared to pQCD calculations. A possible explanation is presented 
in O 
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